Quantum process discrimination with information from environment
Wang Yuan-Mei1, Li Jun-Gang1, †, , Zou Jian1, Xu Bao-Ming2
School of Physics, Beijing Institute of Technology, Beijing 100081, China
School of Physics, Qufu Normal University, Qufu 273165, China

 

† Corresponding author. E-mail: jungl@bit.edu.cn

Project supported by the National Natural Science Foundation of China (Grant Nos. 11274043, 11375025, and 11005008).

Abstract
Abstract

In quantum metrology we usually extract information from the reduced probe system but ignore the information lost inevitably into the environment. However, K. Mølmer [Phys. Rev. Lett. 114, 040401 (2015)] showed that the information lost into the environment has an important effect on improving the successful probability of quantum process discrimination. Here we reconsider the model of a driven atom coupled to an environment and distinguish which of two candidate Hamiltonians governs the dynamics of the whole system. We mainly discuss two measurement methods, one of which obtains only the information from the reduced atom state and the other obtains the information from both the atom and its environment. Interestingly, for the two methods the optimal initial states of the atom, used to improve the successful probability of the process discrimination, are different. By comparing the two methods we find that the partial information from the environment is very useful for the discriminations.

1. Introduction

The overlap between quantum states is a fundamental quantity in quantum mechanics.[16] The overlap between two nonorthogonal states |ψ0〉 and |ψ1〉, O = 〈ψ0|ψ1〉, is useful for a variety of physical phenomena and quantum metrology tasks. For example, it can be used to obtain an upper bound on the probability of successful cloning of the quantum state.[79] Another important application of the overlap is that it can be considered as the essential quantity restricting the effectiveness of discrimination between two states. There are two common approaches for state discrimination.[10,11] The first approach is known as the minimum error state discrimination[1] which gives the minimum error probability, of making an error when two states |ψ0〉 and |ψ1〉 are discriminated, with priori probabilities p0 and p1 = 1 − p0. The second one is the unambiguous state discrimination,[1215] in which the observer is not allowed to give an error instead an inconclusive result. This method can be further optimized by minimizing the probability Pinc of occurrence of the inconclusive result and the minimum probability is So it is obvious that the overlap of two states plays a dominating role in both methods mentioned above.

It is worth stressing that the overlap between two quantum states has an intrinsic property that is preserved when the two states experience the same unitary transformation simultaneously. However, if the states undergo different quantum operations, the overlap between the states will be changed. This property can be used to discriminate quantum operations.[1621] It has been shown that the entanglement of input states[17] or quantum-memory[20] can improve the operation discrimination. Experimental results concerning quantum operations discrimination have been reported.[22,23] These results motivate people to consider the question of whether the overlap can be used to discriminate quantum dynamic processes. The answer is yes. Very recently, reference [24] provided a quantum circuit model to determine which of two candidate Hamiltonians governs the evolution of a quantum system.

We note that realistic quantum systems used in the quantum task are open to essentially uncontrollable environments that act as sources of decoherence and dissipation.[2527] Then the information about an open quantum system will continuously leak out into its environment. According to traditional open quantum system theory,[2527] we always make use of the information that is effectively contained in the system alone and ignore the information that has leaked out into the environment. In some cases, however, it can be possible to extract more information about the state from the environment.[24,28,29] It has been demonstrated that the information obtained by observing an environment that leads to spontaneous decay of a sensing qubit can be used to allow for error correction.[28] Reference [29] modeled noisy channels by coupling the Markov output to “environment” ancillas, and considered the scenario where the environment is monitored to increase the quantum Fisher information of the output. In the quantum process discrimination, it has been shown that additional environment information may supplement what has been obtained from the measurements performed on the reduced system and may therefore improve the successful possibility of distinguishing the processes.[24]

In this paper, we reconsider the process discrimination given in Ref. [24]. In this model the input state of a quantum probe system can be chosen freely and thus we will focus on the question of which initial state of the probe system is the optimal one for the discrimination. For this purpose, we first give a general Hamiltonian of the system driven by a laser field, which can also be used in the non-resonant case. We show that the information lost into the environment has an important effect on improving the successful probability of the discrimination. We also show that when we monitor both the environment and the probe system the excited state is the optimal initial state of the probe system for both resonant and non-resonant cases. However, when we monitor only the probe system, its optimal initial state is always the ground state.

The rest of this paper is organized as follows. In Section 2 we give two measurement methods, one of which obtains only the information from the reduced system state and the other obtains the information from both the system and its environment. In Section 3 we describe in some detail the physical model employed throughout the paper. Section 4 contains our analytical results in the case of resonance and the results of our numerical simulations in the case of non-resonance with our physical explanation about them. Section 5 closes the paper with concluding remarks.

2. General formalism

We suppose that a quantum process is governed by one of two known Hamiltonians H0(t) and H1(t) with prior probabilities p0 and p1 = 1 − p0, respectively. But we do not know which one it is. To discriminate which of the two Hamiltonians is given, we prepare a probe system in a special initial state and then let it access the process. In the dynamic process, the probe system is affected by the environment. Finally, we perform a measurement to determine the process. Here we mainly consider two different measurement methods illustrated in Fig. 1.

Fig. 1. Illustrations of the evolution of a probe system S and the environment E governed by a Hamiltonian H0(t) or H1(t) with a prior probability p0 or p1. In panel (a), we perform an effective measurement on both probe system and environment. While in panel (b) we perform the measurement purely on probe system S.

Method I is shown in Fig. 1(a). In this method we monitor both the environment and the probe system. Suppose that the probe system and its environment are prepared in a state |ψSE(0)〉, and evolves into state or under the Hamiltonian H0(t) or H1(t), then at the final time we can determine which of the candidate Hamiltonians is correct by performing an associate effective measurement to gather the information from both the probe system and the environment as shown in Fig. 1(a). Obviously, the overlap

determines our ability to discriminate the candidate Hamiltonians H0(t) and H1(t). Note that when we calculate the overlap OSE, the states and must be in the same interaction frame.[26] Besides, the overlap can be optimized by choosing an optimal initial state. It is prohibitively complicated to calculate the joint quantum state of the system and environment as it requires one to deal with an infinite number of degrees of freedom. The overlap OSE, however, can be obtained directly without necessarily calculating the whole states. We assume that {|SE〉} is a set of complete bases of the Hilbert space of the probe system and environment. Substituting the complete condition ∑SE |SE〉〈SE| = I into Eq. (1), we can obtain

where

We note that is in the reduced Hilbert space of the probe system and we can obtain the dynamics of the overlap OSE from the time evolution of the . So in what follows, we will give the time evolution of through the formulation of specific measurement schemes in terms of a quantum operation. We assume that the Markov approximation is used, then has a Kraus-like representation[30]

Here, the Kraus operators {Mm} satisfy the normalization condition Suppose

Here, denotes the system Hamiltonian and denote the corresponding Lindblad operators. Substituting Eq. (5) into Eq. (4), expressing and equating term of order dt, we obtain the time evolution equation of

which was first given in Ref. [24] and its alternative derivation can be seen in Refs. [31]–[33]. Equation (6) is structured like the Lindblad master equation, but all operators multiplied with from the left (right) pertain to If we choose we can reduce Eq. (6) to a usual master equation which conserves the trace of the density matrix. Here is not a density operator and the trace of gives the overlap OSE at time t. Then by solving Eq. (6) and taking the trace of we can obtain the time evolution of OSE. Then the distinguishability of the two processes can be evaluated by the minimum error probability

There are two points that should be mentioned: firstly, the two candidate Hamiltonians in Eq. (6) must be in the same interaction frame. Secondly, there is freedom for us to choose an appropriate initial state of the probe system during the process in solving this linear differential equation.

Method II is a traditional one shown in Fig. 1(b). In this method only a measurement on the probe system is allowed at the end of the interaction time t. Then the process discrimination is reduced to distinguish the two density matrices ρ0 and ρ1 evolved by the following Lindblad master equation with the different Hamiltonians,

where ϑ = 0 or 1. The minimal probability of making an assignment error is provided here by Helstrom,[1,26,34] where the summation is over the negative eigenvalues of the operator p1ρ1p0ρ0.

3. Model

Now, we consider a two-level atom driven by a classically single mode laser field. In the Schrödinger picture, the interaction between the atom and the field can be described by the following Hamiltonian in the dipole approximation (ħ = 1):

where |g〉 and |e〉 represent the ground and excited states of the atom, respectively; ωe and ωg are respectively the circular frequencies of the excited and ground states of the atom; ωL and ε0 are respectively the frequency and amplitude of the field; σ+ = |e〉〈g| and σ = |g〉〈e| are the raising and lowering operators, respectively; is the matrix element of the electric dipole moment; φ is the phase of the dipole matrix element. Here, we assume that is a real value and φ = 0. In this model, our quantum process discrimination is essentially to detect which of the fields drives the atom. It is obvious that the field is determined by its frequency ωL, amplitude ε0 and . In general, for convenient calculation, we often transform Hamiltonian Eq. (8) into the interaction frame.

Traditionally, as the first transformation method described in Appendix A, Hamiltonian Eq. (8) in the interaction frame becomes

where δ = ωeωgωL is the detuning of the atomic transition frequency and the driving field frequency. is the Rabi frequency and characterizes the coupling between the atomic dipole and the laser electric field. So in the interaction frame, the field can be described by a pair of δ and Ω. Equation (9) is a time-independent Hamiltonian and convenient for calculation.

Unfortunately, the unitary transformation operator from the Schrödinger picture to this interaction frame depends on the detuning δ (see the Appendix A). Then it will be inconvenient to calculate the overlap OSE from Eq. (9) in the case that there is a nonzero detuning between the atom and the laser field. In order to avoid this problem, we use the second transformation method mentioned in Appendix A, for which the unitary transformation operator does not depend on δ, and obtain another equivalent form of Hamiltonian Eq. (8) in the interaction frame

Here, we note that this Hamiltonian is time-dependent.

Again, we should emphasize that the aim of this paper is to discriminate which of the candidate Hamiltonians governs the quantum process. For notation, we drop subscript I and introduce subscript ϑ = 0 or 1 to represent different Hamiltonians with a prior probability p0 or p1 = 1 − p0. Then equation (10) can be written as

Substitute Eq. (11) into Eq. (6) and choose a single Lindblad operator where κ is the decay rate of the atom excited state decays to its ground state by fluorescence emission, then we will be able to obtain

by which we can obtain and then the overlap OSE. However, according to different laser field frequencies, the quantum process discrimination can be divided into resonance and non-resonant cases and we will discuss these two cases in detail later in this article.

4. Influence of initial state of the atom on the distinguishability

In this section we will discuss the influences of the initial state of the atom on the distinguishability of quantum processes in various cases. We suppose that the two candidate Hamiltonians H0 and H1 have the same prior probability, i.e., p0 = p1 = 1/2. We choose two kinds of initial states of the atom, one kind of initial states are pure states,

and from Eq. (13) we know that θ = 0 or θ = π/2 respectively corresponds to the initial state of the atom being the excited state or ground state. Another kind of initial states are mixed states,

where p is the population of the excited states of the mixed state. It is obvious that the atomic initial state will be the maximally mixed state if p = 1/2. Due to the linearity of Eq. (12), we can deal with the mixed initial state. For this purpose we need to define the overlap of two mixed states as

where

and ε = e, g. represents the overlap between the two states evolved from pure states, either the excited or the ground state.

4.1. Resonance cases

In this subsection we consider the resonance case in which the detuning δ0 = δ1 = 0. Firstly, we consider method I by solving Eq. (12) straightforwardly with and and then evaluate the modulus square of the overlap between states and . The atom is initially prepared in the state Eq. (13) or mixed state Eq. (14). After some calculation we find that the analytic expression of the modulus square of the overlap can be written as

where

and ΔΩ = Ω1Ω0, α = ρeg + ρge with ρi j (i j = eg or ge) being the density matrix element of the initial state of the atom. In the limit ΔΩ → 0.5κ, ε → 0, equation (16) can be replaced by

From Eq. (16), we can find that OSE depends on the initial state of the atom. If the atom is initially prepared in the ground state, excited state or mixed state Eq. (14), then α = 0 and the overlap can be reduced to . We note that is independent of the initial state of the atom and is always real no matter whether ΔΩ > 0.5κ or ΔΩ < 0.5κ. From Eq. (16) we also find that the modulus square of OSE achieves its optimal value if the second term in Eq. (16) vanishes. This can be achieved by initially preparing the atom in the ground state, excited state or mixed state Eq. (14). Then we can come to the conclusion that these initial states are all the optimal initial states of the atom for the Hamiltonian discrimination. While for the initial superposition states (Eq. (13) with θ ≠ 0) the second term on the right-hand side of Eq. (16) will be present then the modulus square of OSE will become larger. When ΔΩ = 0.5κ, we can still obtain the same conclusion that the optimal initial states of the atom are the ground state, excited state or mixed state.

In Fig. 2 we give the time evolutions of the modulus square of the overlap for ΔΩ = 0.3κ (green triangle), ΔΩ = 0.5κ (blue dot), and ΔΩ = 3.0κ (red square). Here we have chosen the excited state as the initial state. While by using the ground state or mixed state mentioned above we can obtain the same results. Figure 2 shows that below the critical value ΔΩ = 0.5κ, the modulus square of the overlap has an exponential decay with time. However, it will vibrate periodically with time for ΔΩ > 0.5κ. From Fig. 2, we also find that the modulus square of the overlap decreases with the increase of ΔΩ. This is because for the resonant case, the difference between the two processes described respectively by H0 and H1 only depends on the Rabi frequency. The larger the ΔΩ, the more easily the H0 and H1 are distinguished from each other. Besides, we find that the modulus square of the overlap always achieves zero in the long time limit, which means that we can distinguish the two candidate Hamiltonians from each other perfectly. Furthermore, from Eq. (16) we can find that the modulus square of the overlap only depends on the value of ΔΩ but is independent of the values of Ω0 and Ω1.

Fig. 2. Time evolutions of the modulus square of the overlap |OSE|2 for ΔΩ = 0.3κ (green triangle), ΔΩ = 0.5κ (critical value, blue dot), and ΔΩ = 3.0κ (red square).

For method II, in which only a measurement on the probe system is allowed, it is not easy for us to obtain the analytic results and thus we give numerical simulations. In Fig. 3 we give the time evolutions of the error probability Ps for different values of Ω0 and Ω1 and different initial states of the atom. The parameters are chosen as follows: in Fig. 3(a) Ω0 = 2.0κ, Ω1 = 5.0κ; in Fig. 3(b) Ω0 = 1.0κ, Ω1 = 4.0κ; in Fig. 3(c) Ω0 = 3.5κ, Ω1 = 4.0κ; in Fig. 3(d) Ω0 = 3.7κ, Ω1 = 4.0κ. We choose the ground (g) state, excited (e) state, superposition (s, θ = π/4) state and maximally mixed (m) state to be the initial states of the atom. The results of numerical calculation show that the error probability Ps depends not only on ΔΩ but also on the concrete values of Ω0 and Ω1, which are illustrated in Figs. 3(a) and 3(b). Besides, we also find that the ground state is the optimal initial state of the atom for large ΔΩ, which is displayed from Figs. 3(b) to 3(d). While when ΔΩ decreases, Ps will increase and the advantage of preparing the atom initially in its ground state will no longer be obvious. We can see that the two curves, which describe the behaviors of Ps for the atom being initially in the maximally mixed state or the superposition state, coincide, which means that the effects that these two states are regarded as the initial states are the same. Another different point is that the error probability obtained by method II does not approach to zero in the long time limit.

Fig. 3. Time evolutions of error probabilities Ps for four initial states of atom in resonance. The four initial states are ground state (g, black square), excited state (e, red dot), superposition state (s, blue triangle) and maximally mixed state (m, green star). The Rabi frequencies are different in the different panels: (a) Ω0 = 2.0κ, Ω1 = 5.0κ, (b) Ω0 = 1.0κ, Ω1 = 4.0κ, (c) Ω0 = 3.5κ, Ω1 = 4.0κ, and (d) Ω0 = 3.7κ, Ω1 = 4.0κ.

These differences between the two methods can be explained as follows. For method I, we consider the information from both the atom and environment while for method II we only consider the information from the reduced state of atom. So the information obtained from the environment is responsible for the improvement of the successful probability of the discrimination.

4.2. Non-resonant cases

In what follows, we will discuss the non-resonant case, i.e., δ ≠ 0. Generally, the quantum process discrimination about our model is to detect which of the laser fields is present and the candidate Hamiltonians As we are unable to give an analytical solution, in this case, we will use a numerical calculation to illustrate our results. We suppose that the laser field is present in only one of the two channels, that is to say, the candidate Hamiltonians and has a given Rabi frequency Ω1 = 4.0κ and adjustable detuning.

Firstly, consider method I. In Fig. 4, we show the modulus squares of the overlap between the two final states for four different initial states of the atom and different values of δ1 = 0.1κ, 0.5κ, 1.5κ, 3.5κ. We find that for small δ1, the dynamic behavior of |OSE|2 for the atom initially in the ground state is similar to that for the atom initially in the excited state or the maximally mixed state. All of them are better than choosing the superposition state as the initial state for most of the time. However, this property mentioned above will not apply to those conditions in which δ1 increases and the advantage of choosing the excited state of the atom as its initial state becomes more obvious. We still note that the maximally mixed state is better than the ground state when δ1 is larger although it is not always the best initial state. Comparing Fig. 4(d) with Figs. 4(b) and 4(c), we find that |OSE|2 will become large when δ1 increases to a certain value. This can be illustrated as follows. The parameter δ1 is the detuning between the atom transition frequency and the frequency of the laser. When the value of the detuning δ1 increases, the effective coupling between the atom and laser decreases, and thus the influence of the driven field will become smaller. This will increase the difficulty in discriminating the process.

Fig. 4. Modulus squares of the overlap for the four initial states of atom when we choose δ0 = 0, Ω0 = 0, Ω1 = 4.0κ and different values of δ1: (a) δ1 = 0.1κ, (b) δ1 = 0.5κ, (c) δ1 = 1.5κ, (d) δ1 = 3.5κ. The initial states are ground (g, black square), excited (e, red dot), superposition (s, blue triangle), and maximally mixed (m, green star) states.

Secondly, we consider method II so as to offer a vivid comparison with our previous method. The results are shown in Fig. 5 which depicts the time evolutions of Ps for different values of δ1. Due to the interaction with the environment, the oscillations are damped over time, and it is intuitively clear that it is no longer optimal to choose the excited state as the initial state of the atom. On the contrary, the ground state is the optimal choice. For large δ1, the minimal value of Ps is increased and the reason for this is the same as that given in the end of the previous paragraph.

Fig. 5. The values of error probability Ps of distinguishing whether a two-level atom is evolved by or initially in the ground (g, black square), excited (e, red dot), superposition (s, blue triangle), and maximally mixed (m, green star) states for Ω1 = 4.0κ and (a) δ1 = 0.1κ, (b) δ1 = 0.5κ, (c) δ1 = 1.5κ, and (d) δ1 = 3.5κ.

Comparing Fig. 4 with Fig. 5, we find that the optimal initial states of the atom are different for the two different methods. The optimal initial state is the ground state for method II, while for method I it should be the excited state. Why is it like this? We can explain it as follows. First we consider the case without energy dissipation, i.e., the decay rate κ = 0. Suppose that the atom is initially in its ground or excited state, following the Schrödinger equation, it will evolve into or after a time t and the analytic expression of the overlap between these two states is

where β = δ1t/2, and The symbol “+” in the above equation corresponds to the atomic initial state that is the ground state while “−” corresponds to the excited state. So it is easy to prove that the distinguishability of the candidate Hamiltonians H0 and H1 are identical no matter whether the atom is initially in the ground state or excited state for the same parameters.

Then a further study is to take the environment into account. In this case, the atom will evolve following the master equation (7). Here we want to reveal the influence of a single quantum channel on the atom state evolution rather than discriminate two quantum processes. So we describe this channel by Hamiltonian HS in this paragraph. In Fig. 6(a), the atom is initially in the ground state. The population of the excited states will be zero if it evolves under HS = 0 and there is no information leaking into the environment. While the population of the excited states will suffer periodical oscillations with damped amplitudes if the atom goes through a channel described by the Hamiltonian HS ≠ 0. At this point, a partial information about the atom will transfer into its environment as the effect of the driven field. Figure 6(b) shows the other cases that the atom is initially in the excited state. We note that the population of the excited states will decay exponentially even when the Hamiltonian HS = 0 and its environment will have some information about the atom. Comparing these two cases, we find that the advantage of the excited state is precisely due to the information extracted from the environment.

Fig. 6. The time evolutions of the population of the excited states of a driven atom coupled to the environment. The initial states of the atom are respectively in (a) ground state and (b) excited state. We consider two cases that the parameters: HS = 0 (orange square), and HS ≠ 0, for δ = 0 (blue dot), δ = 0.5κ (green triangle), δ = 1.5κ (pink circle), δ = 3.5κ (red star), and Ω = 4.0κ.

The above discussion shows that the excited state is always better than the ground state in process discrimination for method I in the non-resonant case. This stimulates us to consider the influence of the population of the excited states of the atomic initial state on the overlap. We compare the ground state with the mixed state as the initial state of the atom for Ω0 = 0, δ0 = 0, δ1 = 1.5κ, and Ω1 = 4.0κ. The results are illustrated in Fig. 7, where p is the initial population of the excited states of the mixed state. We find that as long as the mixed state in Eq. (14) has a nonzero population of the excited states, the modulus square of the overlap compared with the ground state can be reduced. It is because the atom will also leak information into the environment if it initially has a nonzero population of the excited states even though there is no external driven field. It is this piece of information that determines the importance of the population of the excited states in the initial state of the atom. This is a very useful and practical result in an experiment as it is easier for us to prepare such a mixed state than some pure state accurately.

Fig. 7. Modulus squares of the overlap of different initial states, with the initial population of the excited states p = 0 (black), 0.05 (red), 0.1 (blue), and 0.15 (green) (from top to bottom) for Ω0 = 0, δ0 = 0, δ1 = 1.5κ, and Ω1 = 4.0κ.
5. Conclusions

In this work, we reconsider the problem of process discrimination by using the dynamics of the overlap and show that retrieving the information lost in the environment can improve the distinguishability. Specifically, we analyze the optimal initial state of the atom for process discrimination. We mainly discuss two measurement methods. In method I, we monitor the probe and its environment while in method II we monitor only the probe. The difference between these two methods lies in whether the information leaking into the environment is taken into account. For method I we find that in the resonance case, the ground state, excited state or mixed state with zero off-diagonal density matrix element are all the optimal initial states of the atom. However, the optimal state is the excited state in the non-resonance case. Besides we also find that the mixed state mentioned above is better than the ground state and the overlap will decrease as the initial population of the excited states increases. On the contrary, in the case of method II, the optimal initial state of the atom is always the ground state in both the resonance case and the non-resonance case. According to the comparative analyses of the two methods, we can conclude that the information lost into the environment not only contributes to improving the successful probability of process discrimination, but also has an important effect upon the selection of the initial state of the probe system. We hope that our results can stimulate further research on the discrimination of quantum processes.

Reference
1Helstrom C W1976Quantum Detection and Estimation Theory, Mathematics in Science and EngineeringVol. 123New YorkAcademic Press404240–2
2Fox A M2006Quantum Optics: An introductionOxfordOxford University Press35
3Scully M OZubairy M S1997Quantum OpticsCambridgeCambridge University Press55
4Susskind LFriedman A2015Quantum Mechanics: The Theoretical MinimumPenguin BooksUK72
5Barnett S M2009Quantum InformationNew YorkOxford University Press32
6Nielsen M AChuang I L2000Quantum Computation and Quantum InformationNew YorkCambridge University Press409
7Scarani VIblisdir SGisin NAcín A 2005 Rev. Mod. Phys. 77 1225
8Chefles ABarnett S M 1998 J. Phys. A: Math. Gen. 31 10097
9Sekatski PSkotiniotis MDür W 2015 Phys. Rev. 92 022355
10Barnett S MCroke S 2009 Adv. Opt. Photon. 1 238
11Bae JKwek L C 2015 J. Phys. A: Math. Theor. 48 083001
12Ivanovic I D 1987 Phys. Lett. 123 257
13Peres A 1988 Phys. Lett. 128 19
14Chen L YWang H FZhang SYeon K H 2013 Chin. Phys. 22 050306
15Chen L BJin R BLu H 2008 Chin. Phys. 17 778
16Acín A 2001 Phys. Rev. Lett. 87 177901
17D’Ariano G MLo Presti PParis M G A 2001 Phys. Rev. Lett. 87 270404
18Sacchi M F 2005 Phys. Rev. 71 062340
19Zhou X FZhang Y SGuo G C 2007 Phys. Rev. Lett. 99 170401
20Chiribella GD’Ariano G MPerinotti P 2008 Phys. Rev. Lett. 101 180501
21Duan R YFeng YYing M S 2009 Phys. Rev. Lett. 103 210501
22Laing ARudolph TO’Brien J L 2009 Phys. Rev. Lett. 102 160502
23Orieux ASansoni LPersechino MMataloni PRossi MMacchiavello C 2013 Phys. Rev. Lett. 111 220501
24Mølmer K 2015 Phys. Rev. Lett. 114 040401
25Breuer H PPetruccione F2002The Theory of Open Quantum SystemNew YorkOxford University Press
26Wiseman H MMilburn G J2010Quantum Measurement and ControlCambridgeCambridge University Press
27Carmichael H1993An Open System Approach to Quantum Optics, Lecture Notes in PhysicsBerlin/HeidelbergSpringer-Verlag
28Plenio M BHuelga S F 2016 Phys. Rev. 93 032123
29Catana CatalinGutţǎ M 2014 Phys. Rev. 90 012330
30Preskill J1998Lecture Notes for Physics 229: Quantum Information and ComputationCalifornia Institute of Technology
31Gutţǎ M 2011 Phys. Rev. 83 062324
32Gammelmark SMølmer K 2014 Phys. Rev. Lett. 112 170401
33Catana CBouten LGuţă M 2015 J. Phys. A: Math. Theor. 48 365301
34Fuchs C A1996“Distinguishability and Accessible Information in quantum Theory”Ph. D Thesis (The University of New Mexico, Aibuquweque, NM, 1996), arXiv: 9601020